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Three-dimensional gapless semiconductors with quadratic band touching, such as HgTe, a-Sn, 
or Pr 2 lr 207 are believed to display a non-Fermi-liquid ground state due to long-range electron- 
electron interaction. We argue that this state is inherently unstable towards spontaneous formation 
of a (topological) excitonic insulator. The instability can be parameterized by a critical fermion 
number N^- For N < the rotational symmetry is spontaneously broken, the system develops 
a gap in the spectrum, and features a finite nematic order parameter. To leading order in the 
1/A^ expansion and in the static approximation, the analogy with the problem of dynamical mass 
generation in (2-|-l)-dimensional quantum electrodynamics yields W = 16/[37r(7r — 2)]. Taking the 
important dynamical screening effects into account, we hnd that W > 2.6(2) and therefore safely 
above the physical value of = 1. Some experimental consequences of the nematic ground state 
are discussed. 


I. INTRODUCTION 

Three-dimensional (3D) Fermi systems in which va¬ 
lence and conduction bands touch quadratically at the 
Fermi level form the very boundary between two classes 
of materials: Right at the Fermi level their density of 
states vanishes and the systems can hence be under¬ 
stood as limiting cases of semiconductors in which the 
band gap goes to zero. Away from the Fermi level, on 
the other hand, the density of states increases rapidly, 
and the systems may alternatively be regarded as degen¬ 
erate semimetals with the band overlap region shrunk 
to a single point. Prominent examples for such three- 
dimensional gapless semiconductors are given by a-Sn 
and HgTe, which feature band inversion due to spin-orbit 
coupling. The quadratic band touching (QBT) point, lo¬ 
cated at the center F of the Brillouin zone, is protected 
by crystal symmetry, and in the undoped systems the 
Fermi level is right at the touching point [1]. When the 
rotational symmetry is broken, the band degeneracy at 
the F point can be lifted and the spectrum develops a full, 
anisotropic, and topologically nontrivial gap [2]. HgTe, 
for instance, shows the quantum spin Hall effect in quan¬ 
tum well structures [3] and becomes a 3D strong topolog¬ 
ical insulator under uniaxial strain [3]. The pyrochlore 
iridates i?2lr207 (where i? is a rare-earth element) that 
display a metallic paramagnetic state presumably also 
host a 3D QBT point at the Fermi level 0 [3]. This 
scenario has been employed to explain the anomalous 
low-temperature behavior measured in Pr2lr2 07 [7]. 

The dichotomous classification of gapless semiconduc¬ 
tors is reflected in their concomitant peculiar marginal 
screening, which allows fundamentally new types of 
many-body ground states. Theoretical control over the 
situation can be exerted by employing a 1/A^ expansion, 
where N is the number of QBT points at the Fermi level. 
In the limit of large N, Abrikosov and Beneslavskii found 
a scale-invariant semimetallic ground state with anoma¬ 
lous power laws—a 3D non-Fermi-liquid phase [8]. Re¬ 
cently, their analysis has been rediscovered and extended 
in the context of the pyrochlore iridates [7] . 

In this work, we revisit the problem of the ground 


state of 3D Fermi systems with QBT by deriving and 
exploiting the nonperturbative solution of the Dyson- 
Schwinger equations. Although also controlled by the 
small parameter 1/A^, this enables one to address ques¬ 
tions such as spontaneous symmetry breaking and dy¬ 
namical mass generation. We demonstrate that the 
Abrikosov-Beneslavskii large-iV regime exhibits a lower 
bound Nc below which the non-Fermi-liquid state be¬ 
comes unstable and the system features a symmetry¬ 
breaking phase transition. The value for is found to be 
above the physical = I, rendering the 3D systems with 
QBT insulating and nematic at low temperatures. While 
such a possibility has been suggested earlier by consid¬ 
ering the theory in higher spatial dimensions d > 3 [3], 
so far there has not been a definite demonstration of the 
instability within a fully controlled approximation. Con¬ 
ceptually, our scenario is analogous to the well-known 
situation in (2-|-l)-dimensional quantum electrodynamics 
(QED3), which displays chiral symmetry breaking if the 
number of fermions is less than a certain critical fermion 
number [inHil]- The principal difference is, however, 
the lack of relativistic invariance at low energies, which 
makes the effects of dynamical screening both crucial and 
nontrivial to include. We show, nevertheless, that the 
critical number of four-component fermions Nc is, to the 
leading order of the I/7V expansion, bounded from be¬ 
low by 2.6(2), and thus is significantly larger than the 
relevant physical value of iV = 1. 


II. SEMICLASSICAL PICTURE 

The A^-dependent excitonic transition can be under¬ 
stood heuristically within a simple semiclassical picture. 
At the Fermi level, the point of QBT separates the 
filled valence electron band from the unoccupied con¬ 
duction band. Virtual or thermal fluctuations can in¬ 
duce electron-hole pairs, which interact via an attrac¬ 
tive screened Coulomb potential V{r). In the simplest 
approach, the latter is determined by the density of 
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states p(ep) near the Fermi level [15] . 


III. FIELD THEORY FOR 3D QBT SYSTEM 


V{r) 


(Pp e"-P 


The minimal low-energy Hamiltonian relevant for 
semiconductors with diamond or zincblende crystal 
structure is given by the Luttinger Hamiltonian [2Ql [2T] 


where e is the elementary charge and k is the “back¬ 
ground” dielectric constant arising from transitions be¬ 
tween bands away from the Fermi level m- When the 
band dispersion is quadratic, ep oc the density 

of states vanishes linearly as a function of momentum, 
p{ep) oc N\p\, with N the number of QBT points at the 
Fermi level. At large distances, the form of the potential 
thus becomes oc 1/p in Fourier space or, equivalently, 
V{r) oc 1/r^ in real space. The Coulomb interaction 
is marginally screened, in clear contrast to both metallic 
and dielectric screening. Whether or not an electron-hole 
pair can form an excitonic quasiparticle is determined by 
the spectrum of the corresponding s-wave Schrddinger 
equation for the radial “exciton wave function” 'kexc('f)) 

^ exc (r) = 2mE^e^c(r), ( 2 ) 

where m is the effective band mass, a is a dimensionless 
constant, and we have set h = 1. An exciton bound state 
would be given by the solution of Eq. ([^ with energy 
eigenvalue E < 0. The quantum mechanical Hamilto¬ 
nian with 1 /r^ potential is formally invariant under the 
scale transformation r i—>■ /3r, /3 > 0. A discrete bound- 
state spectrum would inevitably break the scale invari¬ 
ance, and one therefore might expect the spectrum to 
possess scattering states only. This naive expectation, 
however, is correct only for weak interaction, when the 
dimensionless parameter a/N < 1/4 [17]. In this case, 
no exciton formation is possible, and the 3D QBT sys¬ 
tem becomes scale invariant at large distances. This is 
the Abrikosov-Beneslavskii non-Fermi-liquid state laiHi. 
For a/N > 1/4, on the other hand, the Hamiltonian 
needs to be regularized at short distances and does ad¬ 
mit bound states, with energy eigenvalues E < 0 that 
depend on the regularization m- This peculiarity of 
the quantum mechanics of the attractive 1 /r^ potential 
can be interpreted as possibly the simplest realization 
of a quantum anomaly m- In our system, screening 
is suppressed at short distances, which naturally regu¬ 
larizes the potential and prevents the electron-hole pair 
from “collapsing” m- For small values of N, exciton 
formation is thus possible and in fact favored energeti¬ 
cally. The critical number of fermions that separates the 
scale-invariant non-Fermi-liquid phase at large N from 
the excitonic insulating phase at small N would in this 
approximation be given by Nc = 4a. 

In what follows we use the field-theoretical machinery 
to demonstrate that the above picture is qualitatively en¬ 
tirely correct. In particular, we will find that the system 
can be described by a differential equation formally iden¬ 
tical to Eq. 1^, with a value for a that leads to A"c > 1- 


I d 

J.2 


dr 


a 1 

]V/2 


Ho(p) 


1 
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(ai -I- |a2) - 


2a2(p- Jf , 


( 3 ) 


with J as the (4x4) j = 3/2 representation of the angular 
momentum algebra, ai the phenomenological Luttinger 
parameters, and where we assumed spherical symmetry. 
The spectrum then is ep = (ai ± 2a2)p^/{2m), and 770 
describes a QBT point if |ai| < 2 |a 2 |- For large N, 
both spherical as well as particle-hole symmetry are in 
fact expected to be emergent at low energies |7], and for 
simplicity we therefore also set ai = 0 in the following. 
While the assumption of spherical symmetry appears to 
be a valid description of HgTe and a-Sn [T], we should 
note that strong correlations in the pyrochlore iridates 
may induce significant anisotropies |22] . For a quanti¬ 
tative description of the latter systems it may therefore 
be necessary to go beyond our simple rotationally and 
particle-hole symmetric model. This is left for future 
work. 

The Hamiltonian can then be written as [3S] 


5 

'Hoi^ = '^da{p)'la, (4) 

a=l 


where da{p) = p^da{'&, p) are the five real spherical har¬ 
monics for the angular momentum of two, viz., di -\-id 2 = 
(•\/3/2) sin^('d)e^®‘^, d^ + id^ = (•\/3/2) sin(2i?)e®‘^, and 
ds = (3cos^ '&—l)/2, with -& and p as the spherical angles 
in momentum space. The Hermitian 4x4 Dirac matri¬ 
ces 7 i ,...,75 satisfy the Clifford algebra { 7 ^, 7 ;,} = 25ab. 
For convenience, in Eq. Q we have also set the remaining 
effective band mass to m/a^ = I- 

We are interested in the effects of the long-range 
Coulomb interaction, which is mediated by a scalar 
field a, and described by the Euclidian bare action 


S' = / 


{dr + ia + Ho) ipi + 


( 5 ) 


ipi and ipi &re four-component fermion fields, and i = 
1,..., N. Upon integrating out a in the functional inte¬ 
gral, the bare l/r density-density interaction is recovered. 

Ho contains the complete set of five anticommuting 
4x4 Dirac matrices; there is no further matrix left 
that would anticommute with all matrices present in the 
Hamiltonian. An isotropic mass gap, that is usually en¬ 
ergetically favored in systems with Dirac fermions [23] 
and in two-dimensional QBT models |25j . is thus im¬ 
possible in the (isotropic) 3D systems with QBT. This 
leaves as the energetically next-best option the full, but 
anisotropic mass gap A oc |5], which can be 
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understood as nematic order parameter |23j . As a gen¬ 
uine many-body phenomenon, dynamical mass genera¬ 
tion is inaccessible to perturbation theory. Similarly to 
QED 3 [TUHT^ and graphene a possible excitonic 

instability can, however, be revealed by a nonperturba- 
tive solution of the gap equation for the fermion Green’s 
function G{uj,'p), 

G~^{uj,p) = iuj + no{iT)+ J 

xG{v,^V{v-uj,q-:^], ( 6 ) 

involving the full vertex function T(u},p',v,q) and 
Coulomb Green’s function V(uj,p). The latter in turn 
is screened by the fermion polarization, 

xG{uj + v,p + ^G{v,^]. (7) 

Similarly, the vertex function fulfills an equation that 
involves the fermion four-point function, and in general 
this leads to an infinite tower of coupled functional inte¬ 
gral equations—the Dyson-Schwinger equations. Within 
the 1/N expansion, the set of equations can be solved 
successively. The wave-function renormalization has the 
expansion Z{uj,p) = 1 -|-C>( 1 /A^) [TO], and to the leading 
order the fermion Green’s function becomes 

= [*u;-hda(^7a + A(a;,^75]"\ (8) 

where we have assumed the effective band mass to be al¬ 
ready at its renormalized value [53] . In Eq. Q and here¬ 
after we adopt the summation convention over repeated 
indices. The vertex renormalization is also suppressed 
by l/N, T{ijj,p\ = 1+0{1/N) [l0|. In QED 3 , higher- 
order calculations give corrections to of no more than 
25% [mill], and it is therefore reasonable to expect the 
relevant physics in the present case also to be captured 
well already by the lowest-order calculation. 

We first discuss the solution in the intermediate range 
of momenta A{uj,p) ^ p^ e^, in which the gap equa¬ 
tion can be linearized. This assumption will be justified 
a posteriori in the limit N ^ Nc- The screened Goulomb 
potential in this regime is, to leading order in 1/N^ 


and interpolating monotonically between these two limits 
for intermediate x. 


IV. STATIC APPROXIMATION 


In the static approximation screened potential V{u},p) 
and gap function A{uj,p) are approximated by their val¬ 
ues at w = 0. Note that in real space then V{r) oc l/r^, 
as anticipated above. If we furthermore assume A(p) = 
A(|pj), the integral over frequency and spherical angles 
in the gap equation Q gives 

NA{p) = ^ r (12) 

37 r( 7 r - 2 ) max(p, q) 


where A « and A « y^A(O) are the UV and IR cutoffs, 
respectively. This integral equation for A{p) is equivalent 
to the differential equation uni 


d / 2 d 

dp \ dp J ' 37r(7r — 2)N 


A{p) = 0, (13) 


supplemented with the boundary conditions 


lim ( A(p) ) = 0 and lim A{p) < 00 . 

p->oo \ dp J p^o 

(14) 


We note that Eq. (13) is formally identical to the ex- 
citon bound-state equation® for E = 0. The very 
same differential equation ([13]) together with the bound¬ 
ary conditions (14) appear also in the leading-order so¬ 
lution of the Dyson-Schwinger equations in QED 3 [TO] . 
quenched QED 4 [50], and graphene [HHZ]. We hence 
simply adopt the results for the present model. For 
N > Nc = 16/[37r(7r —2)] the unique solution to Eqs. (13) 
and (14) is the trivial one, A{p) = 0. For N < Nc, in 


contrast, the equations admit the oscillatory solution 


A(p) oc sin 

Vp 




(15) 


V-^{uj,p} = 4N f 
J u,q 


v{v -I- w) - daiq)daip + q) 


(i/2 -I- q^) [{v -I- a;)2 + {p + (j)^] 


41 > 


(9) 


where we abbreviated ^ = J Carrying out the 

integration gives 


V{uj,p) 



with the dimensionless scaling function E, 


E(x) 


1677 

3 ( 77 - 2 ) 

47rx 


for x ^ 1 , 
for X » 1 , 


( 10 ) 


( 11 ) 


valid for momenta A(0) ^ e^, and with 


A(0)(xexp(^-^=^=+O((iVe/Ar-l)0)^. (16) 

The higher-order terms can be estimated from the numer¬ 
ical solution of the gap equation before linearization [lOj 
as O {{Nc/N - 1)0) =G + 0 {Nc/N - 1) with C « 3.7. 
The transition as a function of N displays an essential 
singularity, similar to the Kosterlitz-Thouless transition, 
and the mass gap is exponentially suppressed in the vicin¬ 
ity of the transition point. This justifies the linearization 
of the gap equation in order to calculate Nc- 
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V. DYNAMICAL SCREENING 


Equations (101 and (111 imply that screening is not as 
efficient at higher frequencies |a;| S> As compared to 
the static approximation, one therefore expects the value 
of Nc to increase when dynamical screening effects are 
taken into account. We demonstrate next that this is in¬ 
deed the case. The linearized gap equation that includes 
the full frequency dependence of the screened potential 
is given by 


NA{uj,^ = 




-Kq4) 


A(z/,q), (17) 


and it defines a Fredholm eigenvalue equation in the 
space of integrable functions with appropriate bound¬ 
ary conditions. The corresponding integral kernel 
can be symmetrized by rescaling the gap function 
A{uj,p)/yJu}‘^ ^ A{uj,p): 


NA{uj,]^ = [k{u;,p;v,q)A{iy,q), (18) 

with k{uj,p]i^,^ = - 

q|-y/(^2 -I- g^)(a;2 +P^)]i and the spectrum of the 
integral operator defined by k is thus real. Each of its 
eigenvectors represents a nontrivial solution to the gap 
equation, and the critical fermion number Nc is hence 
given by the operator’s largest eigenvalue Pmax- 

A lower bound on the value of Nc can be found by 
recognizing that the Rayleigh quotient 

{A\k\A) X.- -A(w,:p)fc(w,:p;i/,r 7 )A(i/,g) 

= - 

is bounded from above by R{k,A) < /imax = Nc- This 
allows a variational approach in which R(k, A) is max¬ 
imized within a suitable set of test functions. As proof 
of concept, we employ this method first for the static 
approximation, where the corresponding integral ker¬ 
nel is given by k^^\p,q) = 4/[37r(7r — 2) max(p, g)]. In 
this case, we can construct an optimal set of variational 
functions from the known solution ( [T^ . We choose 
Aa\p) = sin[alog(^)]/i^/p, with variational parameter 
a and IR cutoff A. For A/A ^ 1 the Rayleigh quotient is 
maximized when a —>■ 0, and it has the asymptotic form 


R(k(°\Ai%,)o, 


16 


37r(7r — 2) 


1 - 


24 


ln(f) 




(19) 

depicted in Fig. ^ In the limit A/A —?► oo the variational 
method approaches the correct static-approximation re¬ 
sult Nc = 16/[37r(7r — 2)], consistent with the fact that 
our ansatz contains the solution (151. 

In order to take the dynamical screening effects 
into account, we use the related ansatz Aa(uj,p) = 



i/Mx) 


FIG. 1. Rayleigh quotient R for a —>■ 0 as a function of 
the cutoff ratio A/A in static approximation (solid line) and 
when dynamical screening effects are taken into account (data 
points). For the latter, the polynomial fit f{x) = fo + fix + 
fsx^ (dashed line) yields Nc > fo = 2.6(2). 


A^^\\p\)/i/^ + p‘^, chosen such that the static solution 
would be recovered if dynamical screening were again ne¬ 
glected in the kernel. To obtain a lower bound on A/, 
we may furthermore approximate the scaling function in¬ 
volved in k{uj,p; u, q) by its lower bound F(x) > Jq)(x) = 
aVT + , with optimized coefficient a ~ 12.17, allow¬ 

ing us to carry out the spherical integrals analytically. 
The result of the remaining numerical integrations over 
frequencies |a;|, \v\ G (A^, A^) and momenta p, g G (A, A) 
is displayed in Fig. as a function of the cutoff ratio 
A/A, for the optimized variational parameter a —f 0. As 
expected, we hnd the Rayleigh quotient, and therefore 
Nc, to be significantly larger than in the static approxi¬ 
mation. The data points are htted well by a third-order 
polynomial inspired from Eq. (191, enabling us to extrap¬ 
olate to A/A —?► oo. This way we obtain a lower bound for 
the critical fermion number with the dynamical screen¬ 
ing effects included: Nc > 2.6(2), and thus significantly 
above N = 1 IHTl. 


VI. EXPERIMENTAL IMPLICATIONS 


The low-temperature ground state is consequently de¬ 
scribed by the mean-field Hamiltonian "Hmf = 7Io + Ao 75 , 
with Aq = A(0,0) > 0. The spectrum of "Hmf is fully 
gapped, with the minimal gap in the Px-Py plane. In 
fact, Hmf describes a QBT system that appears as if un¬ 
der “dynamically generated” uniaxial strain El III 132], 
with the rotational symmetry being spontaneously bro¬ 
ken. Strained HgTe, however, is well-known as a strong 
topological insulator El HI- Without external strain, 
clean HgTe and its analogs hence become at low temper¬ 
atures strong topological Mott insulators E1I331- Equa¬ 
tion (16) allows to estimate the Mott gap: Aq ^ 


e* exp(—L tt/ ^Nc/N — 1 -I- C) > e»/500 when Nc > 2.6, 
with e* the characteristic energy scale for interaction 
effects. For HgTe and a-Sn it is e* ^ 0(1 meV) [T], 
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but larger for Pr 2 lr 207 [5]. Experimentally, the tran¬ 
sition manifests itself, for instance, in the temperature 
dependence of the Hall coefficient Rn, which turns from 
polynomial form i?H c>c above Tc ~ £*/500 A:b to 

an exponential dependence i?H oc exp(Ao/fcB’r) below 
Tc- Similar behavior is expected for the diagonal part of 
the electrical conductivity, which should also inherit the 
anisotropy of the spectrum when T < Tc- 

In uniform magnetic fields h energetics suggests that 
the dynamically induced “strain” aligns parallel to h, and 
is described by 


■Hmf-i-h ='Ho + Ao 75 + h (Jj cos 0 + sin 9) , (20) 


Here, Jz = 17374 -I- § 7172 , and 9 controls the relative 
strength of the symmetry-allowed cubic Zeeman term. 
'Hyni+h leads to a rich phase diagram, previously inves¬ 
tigated in the context of the QBT system under exter¬ 
nal strain [7]. For small 9 and h > he — Ao//iB two 
Weyl nodes with linear dispersion along the Pz axis and 
quadratic dispersion perpendicular to it emerge in the 
spectrum—a pristine double-Weyl semimetal phase |34j . 
By increasing h even further, the Weyl points are shifted 
away from the Fermi level and the system becomes metal¬ 
lic with coexisting Fermi pockets. Across this transition 
the magnetic Griineisen ratio [35] is conjectured to dis¬ 
play a sign change (Tj. When 9 is larger, the transition for 
increasing h is towards a normal metal, before eventually 
the Weyl nodes again emerge. The gapless non-Fermi- 
liquid state, by contrast, turns into a Weyl (semi)metal 
already at infinitesimal magnetic field and should not al¬ 
low any transition at finite h as long as 9 cannot be tuned 
experimentally |7j. 


VII. CONCLUSIONS 

We argued that the 3D gapless semiconductors with 
N points of quadratic band touching at the Fermi level 
exhibit a universal critical fermion number below 
which the Abrikosov-Beneslavskii non-Fermi-liquid state 
becomes unstable. To leading order in the 1/N expansion 
we find = 16/[37r(7r — 2)] in the static approxima¬ 
tion, and even larger when dynamical screening effects 
are taken into account. From the experience gained in 
the related situation in QED 3 mm, we expect that our 
main qualitative result Ac > 1 is robust upon inclusion 
of higher-order corrections in 1/A. At low temperatures, 
clean HgTe and its analogs consequently should suffer a 
phase transition towards a state with spontaneously bro¬ 
ken rotational symmetry, and with a full but anisotropic 
gap in the spectrum—a topological Mott insulator phase. 

We have here limited ourselves to an isotropic model, 
which appears to be a valid low-energy description for 
HgTe and a-Sn, and it leads to the nematic ordering 
as dominant instability. In the pyrochlore iridates that 
exhibit a magnetic transition, however, quantum fluctu¬ 
ations may induce significant anisotropies [22] . In that 
case, a mass gap that anticommutes with the effective 
(anisotropic) Hamiltonian is possible, in contrast to the 
isotropic case, and hence might be favored. It can be 
understood as an order parameter for the antiferromag¬ 
netic all-in/all-out instability [33], and we therefore be¬ 
lieve that a similar scenario as developed in the present 
work might be relevant also for the transitions found ex¬ 
perimentally in the pyrochlores [36] . 
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